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Motivated by recent photoemission experiments on cuprates, the low-lying excitations of a strongly 
correlated superconducting state are studied numerically. It is observed that along the nodal di- 
rection these low-lying one-particle excitations show a linear momentum dependence for a wide 
range of excitation energies and, thus, they do not present a kink-like structure. The nodal Fermi 
velocity vf, as well as other observables, are systematically evaluated directly from the calculated 
dispersions, and they are found to compare well with experiments. It is argued that the parameter 
dependence of vf is quantitatively explained by a simple picture of a renormalized Fermi velocity. 
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Since the discovery of the copper-based high-Tc super- 
conductors 0, there have been extensive studies, both 
experimentally and theoretically, to understand the ori- 
gin of the superconductivity as well as the unusual nor- 
mal state properties • This vast effort has proved that 
among the many experimental techniques angle-resolved 
photoemission spectroscopy (ARPES) is one of the most 
powerful tools, since it can provide important micro- 
scopic information of the electronic structure of these 
materials . One of the recent key findings by ARPES 
experiments concerns the low-lying electronic excitations 
in the (0, 0) to (tt, tt) direction (nodal direction) of the 
Brillouin zone: (i) the low-lying dispersion shows a kink 
at an energy in the range 50-80 meV from the Fermi 
level 0,0, and (ii) the nodal Fermi velocity shows almost 
no doping dependence within the experimental error bars 
(~ 1.8 ± 0.4 eV-A) for a wide range of hole concentra- 
tions X {0 < X ^ 0.2-0.3) 0. To explain the peculiar 
feature (i), the idea that electrons could be strongly cou- 
pled to other degrees of freedom, such as phonons and 
magnetic fluctuations, has been introduced pj- This can 
naturally explain the appearance of a new energy scale. 
There are already theoretical studies in this direction ■ 
However, the kink structure origin is still under debate. 
Regarding feature (ii), there have been no consensus on 
its origin. We believe that before studying more com- 
plicated models it is important to understand to what 
extend a purely electronic model alone can explain these 
features observed experimentally. This is precisely the 
main purpose of this work. 

In this paper, using a variational Monte Carlo (MC) 
method, the low-lying excitations of a strongly correlated 
superconducting state with d-wave pairing symmetry are 
studied. By directly calculating the excitation spectrum, 
it is found that the low energy one-particle excitations 
in the nodal direction show a linear momentum depen- 
dence, instead of a kink-like structure. Our detailed and 
systematic calculations for the nodal Fermi velocity vf, 



as well as a coupling strength A defined below, reveals 
that in spite of not having a kink-structure in £(k) nev- 
ertheless the doping dependence of vp and A are in good 
agreement with experiments. Moreover, it is shown that 
the doping dependence of vp can be understood quanti- 
tatively as a renormalized Fermi velocity. 

As a canonical model for the cuprates, here we consider 
the two dimensional (2D) t-J model on a square lattice 
described by the following Hamiltonian 0; 

H = J ^ (S, • S, - n^n./i) ~ i ^ (4^c,, , -f H.c. 



(1) 



Here is the creation operator of a spin a{=], |) elec- 
tron at site z, and rii — -I- n.i_| and = \<^i a'^apCi^p 
are the electron density and spin operators. («, j) {{hi)') 
runs over the (next) nearest-neighboring sites, and no 
double occupancy is allowed on each site. This model 
has been studied extensively and found to show a d-wave 
superconducting regime in its phase diagram 0, 0| . 

It is well known that a Gutzwiller projected BCS 
wave function with d-wave singlet pairing provides a sat- 
isfactory variational state for the 2D t-J model over 
a wide range of parameters |lll |. Here we use a 
slightly more complex variational wave function 
defined by ^ VnVgVj\QCS) , where |BCS) 

is the BCS ground state wave function, |BCS) = 

Ilk 1 + /kCkjC^ki l'^)' '-k CT Fourier transform of 
cj^g., Vn the projection operator onto the subspace of 
N electrons, Vg = 11^(1 ~ '^j,T'^»,i) the Gutzwiller 

a Jas- 



projection operator, Pj = cxp {^,^ 
trow factor, and /k = Ak/ (Ck — i?k) with Ak 

Avar (cos fc:r - COSfcy), ^k = ^k - /^var , ^k = -2(cOSfc^ 
cos fc. 



4ivar COS kx COS ky , and 5k 



■ V4; ^ ^k 



At \i: 
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Il3| . The variational parameters, all the independent 
pairs of avan ^var, /^var, and t'^^j., are determined 
by minimizing the variational energy: 

for N even [H. In this 

case I^'Im^) is a spin singlet and has a well-defined to- 
tal momentum zero. Hereafter, the wave function with 
the optimized parameters is denoted by and 
= Also the energy unit (t) and the lattice 

constant a are both set to be one. 

A single hole added to is naturally described 

by "'^) = VN-iVGrjll^lBCS), where is the 
creation operator of the standard Bogoliubov quasipar- 
ticle with momentum k and s—l 15]. The state with a 
single added electron, j'i'^^^^-'), is described in a similar 
way replacing Vn-i by Vn+i- Note that these states 
have sharply defined k, total spin 1/2, and ^-component 
of total spin 1/2. The use of these states is partially 
motivated by the diagonalization of H on small clusters 
- effort that indicated that the low-lying single-particle 
excitations are well described by a renormalized Bogoli- 
ubov quasiparticle state [id - and partially on the recent 
proposal of a similar state jl7| . The variational energies 
for these [N ± l)-electron states are denoted by e'^^^\ 
The single-particle excited states dispersion is thereby 
evaluated using £(k)=£;(^) - e'^~'^^ {e'^^^^ - E^^^) 
for the one-electron removal (addition) spectrum ^] . 
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FIG. 1: One-particle dispersion e(k) in the nodal direction, 
k = (fc, k), for the 2D t-J model with J/t=0.3 and t' /t=-Q.2, 
at a;=0.099. (a): Full dispersion for L=162 (triangles) and 242 
(circles). The one-electron removal (addition) spectrum is de- 
noted by open (solid) symbols. The dashed (long dashed) line 
is a fitting curve of branch A (B) for the L=242 data, using 
up to 3rd order polynomials. The estimated Fermi momen- 
tum and energy are indicated by fcp and -Bp, respectively, (b): 
Same as (a) but focusing on the excitations near Ey- In ad- 
dition to the data for L=162 (open triangles) and 242 (open 
squares), results for L = 1250 (solid squares and circles) are 
also plotted. Dotted lines are a guide to the eye. 



In Fig.n](a), as a typical example, the calculated dis- 
persion £(k) in the nodal direction is shown for the 2D 
t-J model with J/t ^ 0.3 and t' /t = —0.2, at hole den- 
sity x=l — ri=0.099. Here n=N/L and L is the total 
size of the system. In the figure, the one-particle removal 
and additional spectra are denoted by open and solid 
marks, respectively. As seen in the figure, these spectra 
are almost symmetric in energy about the center of the 
dispersion. This might be expected because |^(^='=i)) is 
made of a single Bogoliubov mode (with the projections) 
which shows a symmetric spectrum ±f k |l5| . From these 
results several important quantities are evaluated, such 
as the nodal Fermi momentum (kp = (fcp,fcp)), Fermi 
energy (iJp), bandwidth {W), as well as the nodal Fermi 
velocity at kp (up). 

Let us refer to branch A as the main dispersion and 
to branch B as the "shadow" dispersion, since branch 
A (branch B) consists of all the one-electron removal 
(additional) states inside fc™"^, and all the one-electron 
additional (removal) states outside fc™''. Here, fc™'' is 
the nodal Fermi point of |BCS) with the optimized pa- 
rameters. Although this assignment is quite natural, 
we also calculated the quasi-particle weight directly and 
found that branch A has substantially more weight than 
branch B. The energy difference between k = (0, 0) and 
k = (tt, tt) in each branch naturally defines the band- 
width W = £(7r, tt) — £(0, 0). Next, we fit the data in each 
branch using up to third-order polynomials. As shown in 
Fig. ^a), the fitting is highly satisfactory. The intersec- 
tion of the fitting curves provides the Fermi energy Ep 
and the nodal Fermi point kp From these fitting 

curves, we obtain the nodal Fermi velocity v-p at kp. 

The satisfactory fitting of the dispersions £(k) in the 
nodal direction already indicates that £(k) is a smooth 
function of k and, therefore, it suggests that the state 
used here does not have the kink structure observed ex- 
perimentally. To study this in more detail, we calcu- 
lated the nodal dispersion £(k) on a cluster with L=1250, 
where the number of allowed k-points in the nodal direc- 
tion is 50 and, thus, the momentum resolution (5|k| is 
about 0.18/a [l^. The results are presented in Fig.QJb). 
It is fairly clear that the dispersions in both branches 
are almost linear around Ep, and no kink- like structure 
is seen. If there were a kink structure in the dispersion, 
as observed experimentally, it would not be possible to 
fit the data for both one-electron removal and additional 
spectra using the same straight line. Comparing the data 
for L — 1250 and those for smaller systems shown in 
Fig. n^b), it is apparent that the size dependence of the 
dispersion is small, and therefore we can safely estimate 
quantities such as vp using the smaller systems. 

Figure|21summarizes the x dependence of various quan- 
tities for the 2D t-J model with J/t^O.3 and t'/t=-0.2, 
which is a typical parameter set for the cuprates 
Fig. 121 (a) shows that vp is weakly dependent on x up 
to about a;=0. 1-0.2 - with up ~ 0.8-1.0 t - and then 
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increases with further increasing x. li t 500 meV and 
a ~ 4 A are used, the calculated value of this nearly 
x-independent v-p corresponds to approximately 1.6-2.0 
eV-A, which is compatible with ARPES data within the 
experimental error bars [^. The present results are 
also consistent with recent reports by Paramekanti, et 
al. [23, 1^ , including the overall increasing behavior of 
wf as a function of x E^. 
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FIG. 2: (a) Nodal Fermi velocity vy, (b) bandwidth W, (c) 
nodal Fermi momentum |kp| (solid symbols), and (d) renor- 
malized Fermi velocity Wp (solid symbols) (see text) for the 
2D t-J model with J/t=0.3 and t' at different x. 
For comparison, in (c) the free-electrons nodal Fermi points 
(crosses) are also shown. In (d) the Fermi velocity for the free 
electrons Wp (crosses) and up (open symbols) are also plotted. 



In Fig. 121 (b) and (c), the calculated bandwidth W 
and the nodal Fermi momentum k-p are shown, respec- 
tively. Although k-p approximately follows the result for 
free electrons (with t'), W has a stronger x dependence. 
It should be emphasized that this non-trivial x depen- 
dence of W is caused by strong correlations, which are 
imposed in the wave function by the Gutzwiller projec- 
tion Vg- From Fig,|21(b), it is expected that the effective 
mass in the nodal direction monotonically increases with 
decreasing cc, but it is finite even at x ^ 0. 

Now we show that the x dependence of vp is under- 
stood quantitatively by a simple picture of a renormalized 
Fermi velocity. Since kp is similar to the free electrons 
results (with t'), a natural procedure to follow is to cal- 
culate a renormalized Fermi velocity Wp from the value 
of the free electrons Wp at kp. In Fig. 13 (d), Wp = ^v^ is 
plotted, where 7 = W/Wq is a renormalization constant 
and Wo—8t is the free electrons bandwidth. Clearly, Vp 
can now reproduce vp for almost all the doping range 
studied. 

To support this argument, systematic calculations are 
done for various model parameters, and the results are 
shown in Fig. |31 (a)-(d). It is apparent from the fig- 
ures that Vp can indeed explain the x dependence of vp 



quantitatively for a wide hole-doping range. The main 
features of the J and t' dependences are as follows: (i) 
Vp increases with J [Fig. El (a) and (b)], and (ii) the in- 
creasing tendency of vp with x weakens with \t'\ [Fig. |2| 
(b)-(d)]. These dependences can be explained by the 
renormalized velocity picture as well. While kp does 
not depend on J [Fig. ^i)] and, thus, neither does Vp 
[Fig.EIg)], W does depend on J and becomes larger with 
J [Fig.|2Ie)]. Therefore Vp increases with increasing J. In 
contrast, t' does affect the value of Vp, and has a decreas- 
ing trend with increasing x [Fig.|2fg)]. Thus, Vp shows a 
reduced tendency to increase with x. Moreover, this de- 
creasing trend can cancel the increasing behavior of W 
with X and, as a consequence, v^ can present a rather 
weak X dependence over a wide hole-doping range, as in 
Fig. ISJd) , which agrees qualitatively with the calculated 
dependence of vp. For comparison with experiments, we 
also show Ep for various t' in Fig. O (h). The x depen- 
dence of Ep seems to be stronger for large \t'\ than for 
small \t'\, a trend that has been seen in experiments [23l |. 
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FIG. 3: (a)-(d): x dependence of the nodal Fermi velocity up 
(open marks) and the renormalized Fermi velocity Vp (crosses) 
for the 2D t-J model. The parameters used are indicated in 
the figures, (e)-(h): x dependence of W (e), fcp (f), the nodal 
Fermi velocity of the free electrons at kp (g), and -Bp (h) for 
the different parameters indicated in (f). 



Even though our main goal is to study the low-lying 
excitations of the strongly correlated superconducting 
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the bandwidth W decreases with decreasing doping x 
due to correlation, while the bare Fermi velocity instead 
increases and further changes with t'/t. Therefore this 
effect might be less significant from the theoretical point 
of view, contrary to what was previously assumed. Our 
results indicate that more accurate experiments should 
eventually detect a "non-universal" and weakly doping 
dependent Fermi velocity. 



FIG. 4: (a) vf vs. x for the 2D t-J model with J/t=Q.3 and 
t' /t=—0.2 calculated using the Green function MC method 
with fixed node approximation (solid symbols). For compar- 
ison, the variational estimates shown in Fig. |5| (a) are also 
provided (open symbols), (b) x dependence of A for J/t=0.3 
(open marks) and 0.5 (solid marks) with difi^erent t'/t. 
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state, I'S'jj. ), it is interesting to consider the accuracy 
of our estimate of vp for the 2D t-J model. To this end, 
we have carried out a fixed-node approximation Green 
function MC simulation using j'S'j^^^^'') as guiding func- 
tion J^j . As observed in Fig. 0] (a) , the calculated wp 
compare well with the variational estimate . 

Finally, we show in Fig. 0] (b) a coupling strength de- 
fined as X=Vp/vp — 1 which is the first derivative of 
the real part of the single-hole self energy E' at Ep 
(X=—dT,'(u}) /dLu\u!= F;^) if the momentum dependence of 
S' is assumed weak [isl EH . It is interesting to note 
that A seems to be determined largely by J, not by t', 
for a wide range of x, and becomes weakly depending on 
J and t' for x ^ 0.3-0.35. This monotonically decreas- 
ing behavior with x as well as the value of A are in good 
agreement with experimental estimates [^. 

In conclusion, the low-lying one-particle excitations of 
a strongly correlated superconducting state was studied. 
It was found that the dispersion e(k) in the nodal direc- 
tion shows a linear dependence for a wide range of excita- 
tion energies around Ep (~1.0-1.5i) and, thus, does not 
present a kink- like structure . Systematic estimations 
are made for the nodal Fermi velocity vp directly from 
the calculated e(k). The x dependence of vp as well as 
A are in good agreement with experiments. It is shown 
that the model parameter dependence of vp is quantita- 
tively explained by a simple picture of the renormalized 
Fermi velocity. Our results suggest that the kink struc- 
ture observed experimentally is caused by other degrees 
of freedom not included in our study. Although the cal- 
culated e(k) in the nodal direction has no kink, the esti- 
mated Vp as well as the coupling strength A compare well 
with ARPES experiments. These results lead us to the 
speculation that a major part of the low energy physics 
for the cuprates can still be described mainly by a purely 
electronic t-J-like model |22|. The "universal" Fermi ve- 
locity found in experiments 6] turns out to be explained 
here by a rather accidental compensation of two effects: 
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